a^ 

o Directed Percolation with a Wall or Edee 

Q, 

«Q . Per Frojdhf, Martin Howard^, Kent B^kgaard Lauritsen| 

fNORDITA, Blegdamsvej 17, 2100 Copenhagen 0, Denmark 

ICenter for Chaos and Turbulence Studies, The Niels Bohr Institute, Blegdamsvej 17, 

2100 Copenhagen 0, Denmark 






Abstract. We examine the effects of introducing a wall or edge into a directed 
percolation process. Scaling ansatzes are presented for the density and survival 
probability of a cluster in these geometries, and we make the connection to surface 
critical phenomena and field theory. The results of previous numerical work for a 

C^ ■ wall can thus be interpreted in terms of surface exponents satisfying scaling relations 

generalising those for ordinary directed percolation. New exponents for edge directed 
i-rt ', percolation are also introduced. They are calculated in mean-field theory and measured 

C ' numerically in 2 + 1 dimensions. 
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1. Introduction 

The impact of boundaries on critical phenomena has been the focus of much research in 
recent years (for extensive reviews of surface critical phenomena see [^ |^). In the 
presence of a boundary certain (surface) quantities no longer scale as in the bulk, 
but possess different exponents which are dependent on the boundary conditions. In 
the past, most research has focused on the effects of surfaces in equilibrium critical 
phenomena and far less attention has been paid to boundaries in dynamical systems, 
such as in directed percolation (DP). The DP universality class is thought to describe a 
variety of phase transitions from non-trivial active into absorbing states ^ in processes 
such as epidemics, chemical reactions [||, ^, catalysis 0, the contact process 0, and 
certain cellular automata P, Q. Since all of these physical systems contain boundaries, 
an understanding of surface effects is very important. 

The microscopic rules for bulk (bond) DP in d + 1 dimensions are extremely simple: 
any site at time t may make a connection to any of its 2d nearest neighbours at time t + 1 
with growth probability p. Below a threshold, A = p — Pc < 0, such a process always 
dies, whereas for A > there is a finite probability for survival. At the transition point, 
the system is critical and scales anisotropically, i.e. the correlation lengths in time (||) 
and space (_L) scale with different exponents, C,\\ ~ | Aj^'^ii and ^_l ~ | AI^'^-l, respectively. 

Above the upper critical dimension, dc = 4, these exponents can be calculated using 
a simple mean-field theory. However for d < dc, fluctuation effects become important, 
and hence the computation of the exponents becomes a much harder task. The principal 
analytic technique for this calculation employs the equivalence between DP and Reggeon 



field theory [|I0|. Using renormalisation group techniques, the exponents can then be 
computed perturbatively in an e = d^ — d expansion. These analytic techniques are 
supplemented by simulations and series expansions, which mean that, for example, the 



bulk critical exponents for d = 1 are known rather accurately [|lT|- Nevertheless, an 
exact solution for DP remains an open, and extremely important, problem. 

In this paper we will be exclusively interested in the effects of boundaries on DP 
clusters. In order to isolate their effects, it is convenient to consider a semi-infinite 
system, where the cluster grows from a seed close to the surface. Series expansions 



and numerical simulations |T^ in 1+1 dimensions indicate that the presence of the wall 



alters several exponents. In particular, the percolation probability (order parameter), 

Pi(A)~A^\ A>0, (1) 

scales with an exponent jSi rather than the standard exponent (3 (the subscript '1' refers 
to the wall). However, the scaling properties of the correlation lengths (as given by z/y 
and i'±) are not altered. More surprising, however, is the appearance of an apparently 
integer exponent describing the mean lifetime of a finite cluster in the presence of a wall: 



(t)~|A|--. (2) 

Here ti = z/y — jSi = 1.0002 ± 0.0003 in 1 + 1 dimensions and is conjectured to be 
exactly unity |T^. If true, this would be a remarkable result, since none of the other 
exponents for DP are known exactly, and one even lacks evidence for them to be rational 
numbers. Note, however, that this situation is very different from the case of compact 
directed percolation, where no vacancies within a cluster are allowed. This model is 
relatively simple to solve and most of the exponents, including ri, are integers (see M 
and references therein). 

The purpose of the present paper is to analyse the above results in the context of 
surface critical phenomena. We first of all write down scaling ansatzes for the survival 
probability and cluster density which take boundary effects into account. From this we 
are able to derive the behaviour of, for example, the cluster mass in terms of surface 
(and bulk) exponents. We emphasise that the new exponent jSi describes the scaling 
of activity on the wall, thus Pi is a so-called surface exponent [l|, §]. We next consider 
the appropriate field theory for DP in a semi-infinite geometry. This theory was first 



analysed by Janssen et al. [|T5|, where the appropriate surface exponents were computed 
to first order in e = dc — d using renormalisation group techniques. The introduction of 
an inactive wall results in the so-called ordinary surface transition at the bulk critical 
point, for which sites close to the wall are less likely to be active than those in the bulk. 
In this picture, it is clear that only certain exponents are altered: more precisely, the 
boundary introduces one new independent exponent, while all the bulk exponents remain 
unaffected. By placing a finite seed close to the boundary, however, the distribution 
functions of the emerging clusters become sensitive to the new exponent. Unfortunately, 
these field theoretic methods yield little insight into why ti should equal unity in 1 -|- 1 
dimensions. Hence we conclude that the apparently integer value for ri must be a special 
property of DP in 1 + 1 dimensions, inaccessible to perturbative expansions about dc- 

Finally, we extend our analysis by allowing the wall to have an edge with opening 
angle a. This leads to the introduction of new angle-dependent edge exponents which 
govern the properties of clusters started close to the edge. We solve the corresponding 
mean-field theory and also determine the exponents numerically in 2 + 1 dimensions 
using computer simulations. 

2. Wall Analysis 

In this section we shall discuss the exponents associated with the growth of DP clusters in 
the presence of a wall. Some of our analysis will be similar in spirit to that of Grassberger 
T^ , who analysed the case of ordinary percolation in a semi-infinite system. 



First of all, let us examine the effects of introducing a. d — 1 dimensional wall at 
xj_ = [x = (xj|,x_L = 0)] into a DP process. Note that the labels parallel (||) and 
perpendicular (±) refer here to directions relative to the wall (and not relative to the 
time direction). Consider a cluster arising from a single seed located next to the wall at 
t = 0. The probability that an infinite cluster can be grown from this seed is given by 
the percolation probability (|1]) which scales as A'^^ . Furthermore the probability that a 
surface point at a later time belongs to this infinite cluster scales in the same way. Thus 
/?! is an (independent) surface exponent in analogy with surface critical phenomena for 
equilibrium statistical mechanics |jl|, H (more details can be found in section 3, where 
we will also discuss how the surface scaling A^^ crosses over to the bulk scaling A^). 

For a given bulk universality class (such as that of DP), several surface universality 
classes are possible. In our case, the lattice has simply been cut off and hence there will 
be fewer active points close to the surface. This corresponds to the boundary condition 
for the so-called ordinary transition (for which jSi > (3). The survival probability (the 
probability that the cluster is still alive at time t) has the form 

Pi(t,A) = A'^^^i(t/^ll), (3) 

where the scaling function ipi is constant for t ^ ■Cy HI- Furthermore the presence of the 
wall leaves the scaling of the correlation lengths unaltered, and hence exponents such 
as v\_ and z/y are everywhere unaffected. The mean lifetime of finite clusters (^ follows 
from d^) by averaging t with respect to the cluster lifetime distribution —dPi/dt. As a 
result, 

ri = z/||-/?i. (4) 

However, for z/y < /?i, the leading contribution to it) will be a constant, such that the 
above scaling relation breaks down and is replaced by ri = 0. Note, however, if one 
instead considers a space-time geometry where the wall direction departs from the time 
direction, then all the above quantities will, as usual, crossover to bulk scaling (see also 

Hi)- 

For the density pi of a cluster growing from a single seed located next to the wall 
we make the scaling ansatz 

Pi(x,t,A) = Aft+^/i(x/ex,t/e||), (5) 

where the cluster density is defined to be the coarse-grained average density of active 
points. The factor A^^ comes from the probability that an infinite cluster can be grown 
from the seed, whereas A'^ is the probability that the point (x, t) belongs to this infinite 
cluster (see also [|16|)- The shape of the cluster is governed by the scahng function /i. 
In d^) we have assumed that the density is measured at a finite angle d away from the 
wall (where sin ■(9 = x^/x), and suppressed the ■(9- dependence of /i. In contrast, if the 



density is measured along the wall, ^ = 0, then the appropriate ansatz reads 

Pn(a;,t,A)=A2'3i/ii(x/a,t/e||), (6) 

as we pick up a factor A^'^ rather than A'^ for the probability that (x, t) at the wall 
belongs to the infinte cluster. In 1 + 1 dimensions, ■(? has, of course, no meaning. Instead, 
we have a crossover to pu(t, A) = A^^i/ii(t/^||) close to the wall. We also remark that 
for a seed located a (finite) distance away from the wall, the expressions are more 
complicated, although the above scaling forms (||) and (H) are still applicable for large 
times after a crossover from the bulk scaling p(x,t. A) = A'^^f{x/^±, t/^y). 

By integrating the cluster density (||) over space and time, we arrive at the average 
size of finite clusters grown from seeds on the wall, 

{s)r^\Ar\ (7) 

such that 

i^W + du± = Pi + P + Ji. (8) 

Hence, the surface exponent 71 is related to the previously defined exponents via a 
scaling law that naturally generalises the usual d+1 dimensional hyperscaling relation 

z/|| + dz/x = 2/3 + 7. (9) 



It was noted by the authors of [|T^] that relation (|^) is not fulfilled when the exponents 
for the wall geometry are substituted. Within the context of surface critical phenomena 
this "failing" of hyperscaling is perfectly natural, since (^ is only a relation for 
(the unaltered) bulk exponents. The results of previous numerical simulations with 
a wall [p!2| , p!3| are in fact in very good agreement with the modified hyperscaling 
relation (|^). We also note in passing that another generalisation of hyperscaling has 
recently been proposed, although in the rather different context of a (bulk) model 
with multiple absorbing states fl^, where the exponent jS' in the survival probability 
depends continuously on the density of the initial configuration. Such a generalisation 
of hyperscaling might also apply to recent results in [^ for DP with different fractal 
seeds as initial conditions. 

Besides integrating the density (||) , we can also integrate the density on the wall (|^) 
over the d — 1 dimensional wall and time. This integration yields the average (finite) 
cluster size on the wall, 

(Swall) ~ \A\-^'-\ (10) 

where 

z/|, + (rf-l)z/x = 2/5i + 7i,i. (11) 
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However, in higher dimensions {d ^ 2 being a marginal case) this relation is not fulfilled 
as it would predict a negative 71^1. For this case, 71^1 = 0, refiecting a constant 
contribution to (0), cf. the comment after (^. 

The cluster density also contains information on the connectivity correlations, as it 
is proportional to the probability that the seed at the origin is connected to the point 
(x, t). At criticality, we obtain from (^ the power- law decay 

Pi(x,t) = x-(^^+^)/'^^/i(t/x^), (12) 

where z = i'\\/i'± is the dynamical exponent. This is nothing but the critical surface-bulk 
correlation function with pre-factor a;-(^+^i.o)^ which defines the exponent T]ifi describing 
the power-law decay of correlations between the surface and the bulk. Hence, 

/3i + /3 = z^±(rf + r/i,o), (13) 

which generalises the normal DP relation 

2P = u^{d + ri), (14) 

with f] the anomalous dimension in the bulk. Furthermore, by identifying (^) with the 
surface-surface correlation function, it follows that 

2/3i = z^±(rf + r/i,i), (15) 

where rjii is the anomalous surface dimension. As expected, 2?7i = ?7i,i + rj. 

3. Field-Theoretical Analysis 

We now turn to the field theoretic description of DP with a wall and its connections 
with the above scaling picture. The action appropriate for DP with a wall at x_l = is 
given by |T5| 

J = '-'bulk + '-'surface; \^^) 

S^^i^ = Jd'^xjdt (^^[dt-DV'-A]<p + ^u[4><p'-^'<p]y (17) 

'^surface = J d'^-'x J dt A, 0, 0,. (18) 

Here 'S'buik is simply the action from Reggeon field theory [|l^, where is the local 
activity, is the response field, and where we have defined 4>s = 0(x||,x_l = 0,t) and 
(ps = 0(x||,x_L = 0,t). The surface term in S'surface corresponds to the most relevant 
interaction consistent with the symmetries of the problem and which also respects the 
absorbing state condition. Alternatively we can rewrite the action S in the form of a 
Langevin-type equation for the local activity 0(x, t), 

{dt - DV^ - A)0(x, t) + ^M0(x, tf + r7(x, t) = 0, (19) 

(r7(x, t)) = 0, (r/(x, t)r]i^', t')) = w0(x, t)5^(x - x')(5(t - t'), (20) 



where ri{'x, t) is a Gaussian noise term. The multiphcative factor (/)(x, t) in the noise 
correlator reflects the fact that = is the absorbing state. The presence of the wall 
implies the boundary condition at x^ = of Ddxj^(p\g = As4>s- 

For the systematic analysis of DP below the upper critical dimension, the action 
(16)-(18) remains the more useful description. One can show that, in the limit t —^ oo, 
the system reaches a steady-state, where the order parameter (0(x_l,A)) develops a 
profile in the direction away from the wall: 

(0(x^,A)) = A'3v^(xx/a), (21) 

where the exponent (3 describes the density in the bulk, x± ^ ^±, for which the scaling 
function (f is constant. The angular brackets denote averaging with respect to the 
action (16)-(18). Close to the wall, however, the order parameter scales with a different 
exponent than (3. This is in analogy with surface critical phenomena for equilibrium 
statistical mechanics. The new exponent is denoted by j3i, and for A > it governs the 
other limit of the scaling function ip, giving 

(0(x^, A)) ~ A^^xf -^)/^^ XX « e±. (22) 

It is now a standard procedure to derive the form of the correlation functions within 
the field theory. These expressions involve the same /3-exponents as defined above in 
equations ( ^l|) and (^), and are identical to the scaling forms derived earlier in section 
2 (with the exception of an additional inhomogeneous term which is present for the 
surface-surface correlation function 0). This establishes that the /3-exponents defined 
in the field-theory above are indeed the same as the /3-exponents used earlier in the 
scaling theory, which were defined in terms of a percolation probability, as in equation 
(HI). Note that the vanishing of (|2^) in the limit x^ ^ is simply an artifact of the 
continuum analysis (on a lattice the density on the wall simply scales as A'^^). 

Turning now to other aspects of the field theory, it is also straightforward to show (to 
all orders in perturbation theory) that the correlation length exponents are everywhere 



unchanged by the wall — as are all the exponents in the bulk (see 0, |T5|)- Furthermore 
the surface exponent Pi is the only new exponent introduced by the wall. The critical 
exponents can be calculated in a perturbative e expansion around the upper critical 
dimension dc = 4. Hence, quoting from |T5[, we have (identical to the case of DP 
without a boundary) 

(3=l-l + 0{e'), z.|, = l + ^ + 0(62), ^± = l + {^ + 0{e'), (23) 

where e = 4 — d. These exponents are related via hyperscaling (|^) to 7 governing the 
divergence of the bulk susceptibility (average cluster size) and via (|1^ to 77 governing the 



decay of connectivity correlations at criticality. Furthermore an e expansion calculation 



for the surface exponent /3i yields |T5 

A = ^-^ + 0(e2). (24) 

From the field theory of [^, it is not hard to verify that (||) is the appropriate 
generalisation of (||), relating Pi to 

7i = ^ + | + 0(e2), (25) 

which in terms of the field theory describes the divergence of the surface susceptibility 
due to the application of an infinitesimal bulk field. 



The above results are certainly consistent with the numerical work of refs. p2| , [13 
where z/y , i>± were measured in the presence of a wall and found to be unchanged from 
their bulk values. The behaviour oi (3i in (2^) is also in qualitative agreement with 



the available data. Numerically, however, the value of the exponent Ti was found to 
be extremely close to unity in 1 + 1 dimensions. This contrasts with the above series 
results, which give 

n = 4 + ^ + 0(^') (26) 

using the scaling relation (||). Note that the mean- field value of ri appears to be negative. 
This is not in fact the case — following the discussion after (§), we have ri = in high 



enough dimensions. Nevertheless from (|26D we see that the puzzle of why ti seems to 



equal unity in 1 + 1 dimensions cannot be answered by perturbative expansions about 
dc = 4. Therefore this feature would appear to be a special property of DP with a wall 
in 1 + 1 dimensions. This conclusion is certainly in agreement with the 2 + 1 dimensional 



simulations of ref. 13 



4. Edge Analysis 

We next turn to the case of DP in an edge geometry, where the cluster is started on 
an edge. It has been known for some time that the presence of an edge introduces new 
exponents, independent of those associated with the bulk or with a surface (see [1^ 



for a discussion in the context of equilibrium critical phenomena, or |[T^ in the context 
of percolation). However such edge geometries have not yet (to our knowledge) been 
analysed for the case of DP. 

Consider a system, where we allow the wall to have an edge with an angle a at 
3^11 = x± = 0. Hence, the edge can be viewed as the d — 2 dimensional cross section of 
two d — 1 dimensional walls. By placing the seed next to this edge, the surface exponent 
j3i is replaced by the edge exponent /52(a) (where of course /?i = /32(7r)). Following the 
same arguments as before, the survival probability for a cluster starting from the edge 



has the scahng form 

P2{t,A)=A^^^P2{t/^\\), (27) 

where '?/'2 is constant for t ^ ^y. In other words, the percolation probabihty scales as 
P2(A) ~ A^2. Furthermore, we also have the new scaling ansatz for the cluster density 

P2(r,t,A) = A^^+^/2(r/a,t/e||), (28) 

where r is the radial coordinate in a system of spherical polar coordinates centred on 
2^11 = x_L = 0. This ansatz generalises (|^), i.e. it applies for directions away from the 
edge and the walls. By replacing /3 with Pi or P2, we get the corresponding results for 
the density along the wall or the edge, respectively. Moreover, in analogy with (|^) and 
(1) for seeds on a wall, we obtain the average (finite) size (s) ~ |A|~'''^ of clusters grown 
from a seed next to an edge, by integrating (p8|) over space and time. This yields the 
relation 

i/|l +diy^ = (32 + P + 72- (29) 

Similarly, by integrating the corresponding wall density over the d—1 dimensional wall 
and time, we obtain the average size of cluster activity on the wall due to a seed at an 
edge, (s„aii) ~ |A|"^2,i^ ^j^^ 

u\\ + {d - l)u^ = (32 + Pi + 12,1. (30) 



Let us once more remark that scaling relations such as ( p9| ) and ( pO| ) are only valid as 
long as the predicted 7-exponents are non-negative. Our results indicate that 72 should 
be zero for small enough angles a in any dimension, and the same holds for 72,1 also for 
somewhat larger angles. In principle, we can also define an exponent for the average 
cluster size at the edge by (sedgo) ~ jAj''''^'^^ with z/y + {d — 2)z/j_ = 2/52 + 72,2- However, 
after inspecting our numerical results in the next section, we conclude that 72,2 should 
always be zero, with the possible exception of a close to 27r in c? = 2. We also note that 
the wall geometry in section 2 is a special case, such that 72 (vr) = 71 and 72,1(71") = 71 1, 
whereas 72,2 strictly refers to the edge. 

As before, it is also straightforward to identify the various cluster densities with 
correlation functions between different domains p and q. It follows that 

f3p + (3g = u^{d + r]p^g), 2rip^g = {r]p,p + r^,,,), (31) 

with p, g = (bulk), 1 (wall) or 2 (edge). 

We now proceed to calculate the exponents of this geometry in mean-field theory. 



Much of this calculation can be taken over directly from |]T9| where a similar mean- 
field calculation was performed for the case of an Ising model in an edge geometry. 
The appropriate terms in the action (17), (18) yield mean-field equations with some 
resemblance to those of the Ising case. Nevertheless the presence of the time derivative 
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and a different non-linear term leads to some important modifications. However if, at 
the mean-field level, we are interested in calculating either the equal-time two-point 
correlation function or the susceptibility then we can immediately take over the results 



from [jT9|: 



^0,0 = ^ = 0, r/1,1 = 2, 7/2,2 = 27r/a, (32) 

and 

72 = 1 - 7r/2a. (33) 

Furthermore the correlation exponents i>\\ and i'± are again everywhere unchanged by 
the presence of the edge, and retain their usual bulk values. 

However, a calculation of the density exponent /32 requires an analysis of the non- 
linear term. Hence this exponent will differ from that in the edge Ising model. In our 
case, using equation (19), the mean-field local activity for A > must obey the equation 

DVV + A0 - (n/2)02 = 0, (34) 

with the boundary condition that (p -^ 2A/u as r —>■ oo. The solution has the scaling 
form 

0(r) = (A/n)F(r/a,a). (35) 



As r — > the quadratic term in Eq. (|3^ can be neglected, hence we obtain the Ising 
result, with behaving as r^^-^/^. Using i'± = 1/2 we then have oc /\'^+'^/'^°'^ and hence 

^2 = 1 + 71 /2a. (36) 



As a check, we note that this satisfies (^) and (31) at the upper critical dimension 
dc = 4. Similarly, we naively obtain the mean-field values 72,2 = 272,1 = —ir/a < 0, 
which means that 72,2 = 72,1 = 0, as discussed above. Of course we could become more 
sophisticated and use a field theoretic approach to calculate the fluctuation corrections 
to all these mean-field values in an e expansion around dc = 4. Nevertheless below the 
upper critical dimension, we can still expect the mean-field values, and their dependence 
on the angle a, to be qualitatively followed. 

Let us also mention that these ideas can easily be generalised in 3 -|- 1 dimensions 
and higher, where one for example can consider a cluster originating from a seed at the 
cross-section of three walls. The percolation probability and cluster densities will then 
scale with a new corner exponent /^s satisfying scaling relations analogous to ( p9| ) and 
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5. Simulations 

In this section we report the results of simulations of edge DP in 2 + 1 dimensions 
for opening angles of a = 7r/2, 37r/4, vr and 57r/4. We use 2 + 1 dimensional 
bond directed percolation on a bcc lattice where pc = 0.287338(3), and with bulk 
exponents /? = 0.584(5), u^ = 1.295(6), z/_l = 0.734(5), and the dynamic exponent 



z = i^\\/i^± = 1.765(3) ||20| , pTj . In the simulations we start from one seed located on 
the edge (wall for a = vr) and grow the DP cluster. Typically we average over 100,000 
clusters in order to reduce the error bars to a few percent. 
We measure the average position of activity 

{r') = J^^JdV r'p,{r,t,A) = t'^^ h{tA^^^), (37) 

where r is the distance from the seed and the normalisation quantity N{t, A) is the 
mass of the cluster at time t. Thus the average position yields the dynamic exponent 
z = i'\\/i'±- Our results show that z retains its bulk value in agreement with the 
theoretical prediction. Accordingly, we can use the bulk z value in our further analysis 
in order to obtain better estimates for the /32-exponents. 

Next, we measure the critical survival probability which has a power-law behaviour 

P2(t) -r^'/'^ii, (38) 

obtained from Eq. (P7|). The same power law also describes the number of active sites 
(for surviving clusters) on the edge at criticality as a function of time. We also measure 
the probability of having a cluster of mass s which has the critical scaling behaviour 

p{s)^s~-% (39) 

where r^ = 1 + P2/{j^\\ + di'± — P) (cf. ^^)- In addition, we measure the average number 
of active sites at criticality (averaged over all clusters) as function of time by integrating 
(p8|) over space, 

A^(t) -t'^/^-'^/'^ii-^^/'^ii. (40) 

However if we average only over clusters which survive to infinity, then we have instead 

iV.urv(t)~t"/^-^/^ll. (41) 

By measuring the above quantities we can extract various estimates for the ratio 
^2/1^]] which eventually lead to the estimates for P2 listed in table 0. In tables |^ and 
we list our estimates for T2 and 72. These quantities are obtained by measuring the 
average lifetime (t) and average size (s) for finite clusters for different values of A and 
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then obtaining the exponents by carrying out a power-law fit. We observe that the 
results in 2 + 1 dimensions qualitatively show the behaviour expected from the mean- 
field predictions. With one exception, we confirm that the scaling relation 

r2 = v\\ - P2 (42) 

[cf. the analogous expression (D for a wall] and hyperscaling (^) are both fulfilled when 



error-bars are taken into account. This exception occurs for the smallest angle where 
the relation (^) is not fulfilled. This is because as soon as /32 becomes larger than z/y (as 
is the case for a = vr/2) the above relation breaks down, and instead the mean cluster 
lifetime becomes constant (i.e. T2 = 0), cf. the comment after (^). Using our results for 
/32 we find that in 2 -|- 1 dimensions T2 will reach its mean-field value of zero for an angle 
in between 7r/2 and 3tt/A. When T2 approaches zero the correction to scaling terms in 
the expression for (t) will affect the scaling making it difficult to obtain precise values. 

6. Conclusions 

We have analysed the impact of a wall or edge on a directed percolation process in terms 
of surface critical phenomena. The presence of an inactive wall results in an ordinary 
phase transition between active and inactive surface states at the bulk critical point. 
A description of this transition requires the introduction of one further independent 
exponent in addition to those present in the bulk. We have formulated a scaling 
ansatz for clusters growing near the surface, which incorporates the surface effects and 
explains how the exponent for the survival probability is altered. We have also used 
the connection between DP and Reggeon field theory to justify our scaling ansatzes and 
to examine the nature of the surface exponents. It turns out that the conjecture for 
Ti = 1 in 1 + 1 dimensions cannot be explained within the e expansion. We also remark 
that it would be possible to examine other surface universality classes for evidence of 
rational exponents, particularly at the special transition. This transition occurs when 
the surface bond probabilities are enhanced such that not only the bulk, but also the 
surface, is at criticality. We note that the transition at this (multicritical) point requires 
the introduction of two new independent exponents. Lastly, we have for the first time 
analysed edge exponents in DP for edges with variable opening angles. We have derived 
the mean-field exponents and computed numerical values from computer simulations in 
2 + 1 dimensions. 
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Table 1. Estimates for the (32 exponents for 2 + 1 dimensional edge DP together 
with the mean-field values. The bulk and 1 + 1 dimensional wall estimates [n3, n3| are 
listed for reference. The mean-field value /32 is obtained from Eq. (pq). Recall that 
/32(^) = /3i. 



Angle (a) 


7r/2 


37r/4 


71 


57r/4 


bulk 


pr^ 






0.7338 ±0.0001 




0.2765 ±0.0001 


pr^ 


1.6 ±0.1 


1.23 ±0.07 


1.07 ±0.05 


0.98 ±0.05 


0.584 ±0.005 


^(MF) 


2 


5/3 


3/2 


7/5 


1 



14 



Table 2. Estimates for the T2 exponents for 2 + 1 dimensional edge DP. The bulk and 
1 + 1 dimensional wall estimates [O, 0| are listed for reference. Note that T2 =0. 
Recall that T2{tt) — ti. 



Angle (a) 


7r/2 


37r/4 


vr 


57r/4 


bulk 


^1 






1.0002 ±0.0003 




1.4573 ±0.0002 


^(2+1) 


0.1 ±0.05 


0.20 ±0.05 


0.26 ±0.02 


0.38 ±0.04 


0.711 ±0.007 



Table 3. Estimates for the 72 exponents for 2 + 1 dimensional edge DP together 

with the mean-field values. The bulk and 1 + 1 dimensional wall estimates fn3, Ol are 

rMFi I I — — 

listed for reference. The mean- field value 72 is obtained from Eq. (p3|). Recall that 

72(71") = 71. 



Angle (a) 


7r/2 


37r/4 


TT 


57r/4 


bulk 


,.'.., 






1.8207 ±0.0004 




2.2777 ±0.0001 


7^" 


0.7 ±0.1 


1.0 ±0.1 


1.05 ±0.02 


1.20 ±0.05 


1.592 ±0.009 


,r> 





1/3 


1/2 


3/5 


1 



